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Abstract 

Magnetization measurements in the low field region have been carefully performed on a well- 
shaped cylindrical and an ellipsoidal sample of superconductor MgCNi%. Data from both samples 
show almost the same results. The lower critical field H c \ and the London penetration depth 
A are thus derived. It is found that the result of normalized superfluid density A 2 (0)/A 2 (T) of 
MgCNi$ can be well described by BCS prediction with the expectation for an isotropic s-wave 
sup er conductivity. 
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I. INTRODUCTION 



The pairing symmetry is very essential for uncovering the mechanism both for con- 
ventional and high-T c superconductivity. The recently discovered intermetallic perovskite 
MgCNi 3 superconductor— is regarded as a bridge between conventional superconductors and 
high-T c cuprates, and the issue concerning its symmetry of order parameter has attracted 
considerable attention. However, pairing symmetry about MgCNi 3 remains highly con- 
troversial in reported literatures. NMR— , specific heat£, scanning tunneling measurement^ 
and point contact tunneling spectra^ favor the s-wave pairing in MgCNi 3 . On the other 
hand, the earlier theoretical calculation^, the tunneling spectral and the penetration depth 
measurement 8 support non-s-wave superconductivity. Recently a two-band s-wave model 
has been proposed by Wdlte et al£ who try to explain the complex behavior observed in 
MgCNi 3 . 

In this paper, we derive the thermodynamic parameters H c i and A of two MgCNi 3 
samples by careful magnetization measurement. It is found that the normalized superfluid 
density, A 2 (0)/A 2 (T), can be described by BCS prediction for a s-wave pairing symmetry. 
Therefore, our magnetization data support the conventional single band s-wave supercon- 
ductivity in MgCNi%. 

This paper is organized as follows: The samples and experimental details are presented 
in section II. The data and discussions are given in section III. And section IV gives the 
summary. 

II. SAMPLES AND EXPERIMENTAL DETAILS 

The polycrystalline MgCNi 3 sample investigated here has been prepared by powder 
metallurgy method, and the details of preparation can be found elsewhere^. The supercon- 
ducting transition temperature is 6.9 K measured by both magnetization [ ac susceptibility 
( / = 333 Hz, H ac = 1 Oe ) and dc diamagnetization shown in Fig. 1(a) ] and resistivity 
measurement. The M(T) curves show a sharp transition with the transition width less than 
0.5 K. The x-ray diffraction ( XRD ) analysis presented in Fig. 1 (b) shows that all diffrac- 
tion peaks are from the MgCNi 3 phase, which indicates that the sample is nearly of single 
phase. 
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FIG. 1: (a) The ac susceptibility curves of MgCNi^ measured by MagLab with ac field 1 Oe and 
frequency 333 Hz. The inset gives the ZFC and FC dc diamagnetization at 20 Oe measured by 
SQUID, (b) XRD patterns of the MgCNi 3 . 

In order to minimize the demagnetization factor, one sample ( denoted as S-c ) has 
been carefully cut and ground to a cylinder with a diameter of 1.1 mm and length of 7.0 
mm. The demagnetization factor in this situation is almost negligible since the field has 
been applied along the axis of the cylinder. Another sample ( denoted as S-e ) has been 
polished to an ellipsoid with semi-major axis a=3.74 mm and semi-minor axis b=1.5 mm. 
The demagnetization factor for the ellipsoidal sample is n = (1 — -^)(1 — ^ • In j^-) ~ 0.136, 
with e = \Jl — b 2 /a 2 . The magnetic fields have been applied parallel to the longitudinal 
axis of the samples. 

The magnetic measurements are mainly carried out on an Oxford cryogenic MagLab 
system ( MagLabl2Exa, with temperature down to 1.5 K ) and checked by a quantum 
design superconducting interference device ( SQUID, MPMS 5.5 T ). After zero-field cooled 
( ZFC ) from 25 K to a desired temperature, the magnetization curve M(H) is measured 
with the applied magnetic field swept slowly up to 1000 Oe ( ^> H cl ). It is important to 
note that the magnet has been degaussed at T= 25 K in order to eliminate the remanent 
field before each measurement. It is essential to do degaussing since otherwise even 5 Oe 
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residual field may cause significant effect on the result of magnetization. 

III. EXPERIMENTAL DATA AND DISCUSSIONS 

In this section, the processes to obtain the superconducting parameters by magnetization 
measurement have been reported in detail for two MgCNi^ samples, one is a cylinder and 
another is an ellipsoid. 

A. The cylindrical sample ( S-c ) 

The curves of dc magnetization are shown in Fig. 2. The temperature varies between 
1.57 K and 6.88 K with steps 0.2 K ( some 0.4 K ). All curves show clearly the common 
linear dependence of the magnetization on field caused by Meissner effect at low fields, and 
this extrapolated common line is the so-called "Meissner line"( ML ). The optimal ML ( 
solid line in Fig. 2 ) is achieved by doing linear fit M(H) of the lowest temperature ( 1.57 
K ) at low fields, which represents the magnetization curve of Meissner state. The value of 
H cl is determined by examining the point of departure from linearity on the initial slope 
of the magnetization curve ( ML ) with a certain criterion. The results of subtracting this 
ML from magnetization curves are plotted in Fig. 3 and the AM between 7.0 x 10~ 4 and 
1.4 x 10~ 3 emu are shown in the inset with an enlarged view. All curves show a fast drop 
to the resolution of device when the real H cl is approached, so the value of H cl is easily 
obtained by choosing a proper criterion of AM. The H c i(T) acquired by using criteria of 
AM = 7.0 x 10~ 4 and 1.1 x 10~ 3 emu are shown in Fig. 4. Then the penetration depth 
A(T) can be achieved from H cl (T) by 

and they are displayed in the inset of Fig. 4. Here $ = hc/2e k 2 x 10~ 7 G ■ cm 2 is the 
flux quantum, and k, is the Ginzburg-Landau parameter. We take k as constant since it is 
a weakly temperature dependent parameter. 

The values of nominal H c i and A seem to be criterion dependent in this method, however 
temperature dependence of H c i and A(T) are found to be weekly criterion dependent if the 
data is normalized by the zero temperature values ( see A 2 (0)/A 2 (T) in Fig. 5 ) . In addition, 



4 



200 400 600 800 1000 

H(Oe) 



FIG. 2: The magnetization curves of M(H) for MgCNi^. The solid line is the "Meissner line" de- 
fined in the text. The temperature is varied from 1.57 K to 6.88 K from bottom to top, with step 
about 0.2 K ( some 0.4 K ). It is found that the initial slope of all the curves is the same. 
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FIG. 3: The difference between the "Meissner line" and the M(H) curves. AM is shown in 
logarithmic scale. The inset shows the enlarged AM ( between 7.0 x 10~ 4 and 1.4 x 10~ 3 emu ). 

AM drops sharply with decreasing magnetic field, the use of lower AM value in our criterion 
will not result in a much different H cl (T) curve. If not specially mentioned, the discussion 
is based on the data using the criterion of 7.0 x 1CT 4 emu hereinafter. At the temperatures 
below 2.8 K ( < 0.4 T c ), the values of H c i(T) and A(T) are almost constant despite the 
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FIG. 4: The temperature dependence of the nominal H c \ and A ( inset ). H C \{T) is obtained by 
using criterion of AM = 7.0 x 10 -4 ( open circles ) and 1.1 x 10~ 3 emu ( solid circles ), respectively. 
Error bars are given for determining the nominal H c \ and A. Lines are guides to the eyes. 
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FIG. 5: The temperature dependence of A 2 (0)/A 2 (T) for MgCNi 3 and MgB 2 . The values of 
A 2 (0)/A 2 (T) are obtained by H cl (T) using criteria of AM equal 7.0 x 10 ( open circles ) and 
1.1 x 10~ 3 emu ( solid circles ) for MgCNi 3 and 1.0 x 10~ 4 emu for Mg£>2( solid). The prediction 
by two-fluid model ( dotted ) and BCS s-wave( dashed ) are also shown. The data of open and solid 
triangles are the experimental measurement by microwave resonator method and radio frequency 
technique for MgB2, respectively ( from Eef. and Ref.|17j| ). 

lack of the data below 1.5 K. This may imply the conventional s-wave nature in MgCNi 3 , 
because the finite energy gap manifests itself with an exponentially activated temperature 
dependence of thermodynamic parameters. This can be further confirmed in the following 
discussion on superfluid density. Worthy of noting is that, it is very difficult to distinguish 
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a slight difference of H& ( or 1/A 2 ) in low temperature region between different pairing 
symmetries, for example for an ideal s-wave, an exponential dependence is anticipated, for a 
dirty d-wave, a quadratic form p s (T) = p s (0) — <xT 2 is expected. Here we use an alternative 
way, i.e, to fit the data in intermediate and high temperature region to extract useful message 
for pairing symmetry. 

As we know, the total superfluid density p s is proportional to A~ 2 (T), and A 2 (0)/A 2 (T) 
represents the normalized superfluid density. In Fig. 5, we display the temperature de- 
pendence of A 2 (0)/A 2 (T) of MgCNi 3 with A(0) as a fit parameter. The predictions of BCS 
s-wave( dashed ) and two- fluid model ( dotted ) are also shown. According to the BCS theory 
for clean superconductor a 11 ^ 2 , the normalized superfluid density A 2 (0)/A 2 (T) is expressed 



where A(T) is the BCS superconducting energy gap, f(E) = l/[exp(—E/k B T) + 1] is the 
Fermi distribution function, and D(E) = E/(E 2 — A 2 (T)) 1 / 2 is the quasiparticle density of 
states. The most appropriate superconducting gap A(0) = 1.86A;&T C is chosen in our BCS 
calculation with T c = 6.9K, and this value is reasonable for MgCNi 3 because the generally 
reported results are larger than the conventional BCS value (1. 76 ki,T c ). It is found that 
A 2 (0)/A 2 (T) of MgCNis can be well described by the s-wave BCS theory with a single gap, 
but the two-fluid model shows a substantial deviation. This suggests the s-wave nature of 
superconductivity in MgCNi%, which is consistent with our previous conclusion reached by 
point-contact-tunneling^. Later on we will show that our results are not compatible with 
any other pairing symmetry with nodes on the gap function which normally contributes a 
power law dependence to the temperature dependence p s . 

For the sake of comparison, the temperature dependence of the normalized superfluid 
density in MgB 2 obtained by exactly the same magnetization method^ is also shown in 
Fig. 5, with A(0) as a fit parameter. Clearly the data can not be understood in isotropic 
s-wave BCS theory or two-fluid model because of the two-gap characteristic of MgB' 2 l4]l5 . 
The data obtained from this simple magnetization method on MgB 2 was found to be close 
to that determined by more elegant microwave method 16 and radio frequency technique^, 
which can be seen from Fig. 5. This indicates that the same magnetization method used in 
MgCNi 3 to get H cl and A is reliable, and the corresponding results are plausible. 

One may argue whether the H c i obtained here is the lower critical field of grains because 
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of the polycrystalline nature of our sample. In our magnetization experiment, the nominal 
H c i(0) of our sample is about 145.1 Oe. Combined with H C 2(0) ( 1.18 x 10 5 Oe ) determined 
from our previous measurement of specific-heat^, we can reach that the value of k is 39 and 
H c (0) equals to 2165 Oe by Eqs. (3, 4). The value of coherence length is 5.3 nm obtained 
by Eq. (5). 

H c1 {T) = ±=H c {T)Uuk (3) 
H c2 (T) = V2H c (T)k (4) 



£(0) = ^$ /27Ti/ c2 (0) (5) 

And the value of A(0) is about 200.1 nm. All these values of parameters are in the range of 
the reported results of MgCNi 3 by other techniques( see collected parameters in Ref.|9j ). 
This manifests that H c \ measured here reflects the bulk property. In addition, the value of 
£ for MgCNi 3 is quite large, so that the influence of the grain boundary is weak. 

Another argument is that the nominal H c % relation obtained in our experiment may not 
reflect the true H c \ but the flux entry field because of the Bean-Livingston surface barrier and 
effects of sample corners geometrical barriers. However, we would argue that the influence of 
surface barrier is not important to our cylindrical sample, since the magnetization hysteresis 
loops are very symmetric in the temperature and filed regimes we measured. In order 
to further verify the validity of this method to obtain H&, we have repeated the same 
measurement for an ellipsoidal sample. The data and the discussion are presented below. 



B. The ellipsoidal sample ( S-e ) 

The curves of dc magnetization are shown in Fig. 6 and the temperature varies between 
1.59 K and 6.90 K with steps about 0.1 K . The optimal "Meissner line"( solid line in Fig. 6 
) has been determined in the same way as for the cylindrical sample. Subtracting this ML 
from the magnetization data yields the AM curves plotted in Fig. 7. The H cl (T) acquired 
by using criteria of AM = 3.2 x 10~ 4 and 1.0 x 10 -3 emu are shown in Fig. 8, and the 
demagnetization factor n ( w 0.136 ) has been taken into account. Then the penetration 
depth A(T) can be achieved from Eq(l) and they are displayed in the inset of Fig. 8. The 
normalized temperature dependence of A 2 (0)/A 2 (T) of MgCNi 3 is shown in Fig. 9. One 
can clearly see that the data from the ellipsoidal sample is almost identical to that for the 
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FIG. 6: The magnetization curves of M(H) for the ellipsoidal MgCNi^ sample. The solid line is 
the "Meissner line". The temperature is varied from 1.59 K to 6.90 K ( from bottom to top ), with 
step about 0.1 K . 



FIG. 7: The difference between the "Meissner line" and the M(H) curves. AM is shown in 
logarithmic scale. 

cylindrical sample, showing a trivial influence of either the geometrical or surface barrier in 
our present samples. 

In addition, we have calculated the superfluid density assuming a nodal gap with d-wave 
symmetry. Under the frame of the BCS theory, if the gap has a d-wave-like node, the 
normalized superfluid density A 2 (0)/A 2 (T) is written as 



with f(E) = l/[exp(-E/k B T) + 1], A(T, 6) = A (T) • cos 29 and D(E) = 
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FIG. 8: The temperature dependence of the nominal H c \ and A (inset). Lines are guides to the 
eyes. 
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FIG. 9: The temperature dependence of A 2 (0)/A 2 (T) for ellipsoidal MgCNi% sample. The values 
of A 2 (0)/A 2 (T) are obtained by H c i(T) using criteria of AM equal 3.2 x 10~ 4 ( open circles ) and 
1.0 x 10 -3 emu ( solid circles ). The prediction by two-fluid model ( dotted ), BCS s-wave ( solid 
) and pure d-wave ( dashed ) are also shown. 

E/(E 2 — A 2 (T, 6)) 1 / 2 . The calculated d-wave results are shown in Fig. 9 with dashed line. 
For p x -wave symmetry with A(T, 6) = A (T) • sin#, we found that the calculated tempera- 
ture dependence of A 2 (0)/A 2 (T) is close to that of d-wave, and far from our experiment data. 
The predictions of s-wave BCS ( solid ) and two-fluid model ( dotted ) are also shown in Fig. 
9. It is found that our data can only be well described by the s-wave model. Together with 
the results for the cylindrical sample, we conclude that MgCNi 3 is most likely an isotropic 
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s-wave superconductor. 



IV. SUMMARY 

To summarize, we have measured the M-H curves of two MgCNi^ samples with cylin- 
drical and ellipsoidal shapes and obtained their lower critical field H c i(T) and A(T). The 
temperature dependence of normalized superfluid density is consistent with the s-wave BCS 
theory. All these indicate that MgCNi^ may possess an isotropic s-wave gap, which is in 
sharp contrast to MgB 2 . 

Note added: The recent report of carbon isotope effect in MgCNi3 by T. Klimczuk and 
R.J. Cava indicates that carbon-based phonons play an essential role in the superconducting 
mechanism^. 
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